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Abstract 

We extend a non local and non covariant version of the Thirring model 
in order to describe a many-body system with spin-flipping interactions 
By introducing a model with two fermion species we are able to avoid the 
use of non abelian bosonization which is needed in a previous approach. 
We obtain a bosonized expression for the partition function, describing 
the dynamics of the collective modes of this system. By using the self- 
consistent harmonic approximation we found a formula for the gap of the 
spin-charge excitations as functional of arbitrary electron-electron poten- 
tials. 
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1 Introduction 



In recent years there has been a surge of activity in the study of low- 
dimensional field theories. In particular, research on the one-dimensional (Id) 
fermionic gas has attracted a lot of attention. An interesting aspect of this 
system is the possibility of having a deviation from the usual Fermi-liquid be- 
havior. This phenomenon is known as Luttingcr-liquid behavior, characterized 
by spin-charge separation and by non universal (interaction dependent) power- 
law correlation functions jD. The simplest theoretical framework that presents 
this feature is the Tomonaga-Luttinger (TLM) model a many-body system 
of right and left-moving particles interacting through their charge densities. In 
recent papers || Q] || an alternative, field theoretical approach was developed 
to consider this problem. In these works a non local and non covariant version 
of the Thirring model was introduced, in which the fermionic densities and cur- 
rents are coupled through bilocal, distance-dependent potentials. The resulting 
non local Thirring model (NLTM) contains the TLM as a particular case. Al- 
though it provides an elegant framework to analyze the Id many-body problem, 
as it stands it is not very practical to consider magnetic properties of Luttingcr 
liquids. This is so because in the original formulation of NLTM ||, the descrip- 
tion of spin-flipping (SF) processes, based on non abelian bosonization, is quite 
involved. In this paper we show how to circumvent this problem by introducing 
a Thirring-like model with two fermion species. To this end we generalize the 
two-fermion model of Ref. Q , originally built as a local and covariant theory, 
to the case in which the interactions between fermionic currents are mediated 
by bilocal functions. From the condensed matter standpoint, our work can also 
be viewed as an extension of the work of Grinstein, Minnhagen and Rosengren 
where a simplified version of the SF problem (contained in our model) was 
first considered. Our procedure allows to get a bosonized vacuum to vacuum 
functional by carefully extending previous formulations of massive § [pj and 
massive- like @ models. The paper is organized as follows. In Section 2 we 
introduce the above mentioned non local two-fermion model and explain its re- 
lationship to previous non Abelian description. In Section 3 we establish an 
equivalence between the initial fermionic partition function and the one cor- 
responding to a two-boson non local extension of the sine-Gordon model. In 
Section 4 we employ a non local version of the self-consistent harmonic ap- 
proximation (SCHA) jll| in order to analyze the spectrum of the spin-density 
modes. In particular we give, within this approximation, a closed formula for 
the gap as functional of the non-contact forward-scattering (no SF) potentials. 
Finally, in Section 5, we summarize the main points of our investigation and 
gather our conclusions. In the Appendix we show how to compute the gap of 
the spin-charge sector by means of the harmonic approximation. 
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2 The model and its relationship to a previous 
non Abelian description 

In this Section we introduce a non local version of the Thirring model which 
incorporates electronic spin by considering two fermion species, where each 
species represents a spin state. Our initial (Euclidean) action is 

£x^m a ~ 9 -^ J <Pxd 2 y] «{x)V${x--y)]l{y) 

-g s J dPxiPy&ix^&WUwix-v)* 2 ^* 1 ^) (1) 
where a, b — 1, 2, with 1 ="[", 2 =|, the currents are the usual fermion currents 

jl = jl = * 2 7^ 2 , (2) 

and the matrices V, a \ have the form 



(3) 



where and V? are coupling potentials related to Solyom's .92 and 54 
forward scattering couplings. These functions are related to current-current 
interactions which do not flip the spins. The functions U^(x — y) are the 
couplings governing SF processes. We have kept the constants g and gs in 
order to facilitate comparison with local versions of this model. For instance, 
the case gs — and Vffi) — Vffi) — 5 2 (x — y) corresponds to two decoupled, 
usual Thirring models. The action (Q) has manifest t/(l) chiral symmetry, i.e., 
it is invariant under the transformation ^> a — > e* 7B6, \l/ a , ^ a — > $ a e I75e . Fermion 
spin conservation is also preserved in this theory. 

The theory defined above is similar to the one described by ZJ ||. There 
are, however, two important differences with that model. First of all, our model 
takes into account the possible long-range nature of the potentials, whereas ZJ's 
model is local. On the other hand ZJ's model does not include c/4-like terms 
associated to scattering processes involving just one fcrmionic (left or right) 
branch, neither for SF nor for ordinary diagrams. 

Concerning the relationship of our action to previous condensed-matter in- 
spired models, we should mention the pioneering works of Luther and Emery 
[|l3|, and Grinstein, Minnhagen and Rosengren M. The first authors introduced 
the so called backscattering model. Although this system has in principle no 
SF, in the local limit their backscattering diagrams coincide with the ones of 
spin-changing processes. Grinstein et al. included from the beginning SF in- 
teractions taking into account a Coulombian potential. Thus, their model is 
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non local, but they took the same potential for all kinds of diagrams (SF and 
ordinary). Besides, in order to establish the equivalence between their theory 
and a Coulomb gas system the authors considered again a local limit. These 
authors did not include (74-like terms either. 

Let us now show that the action ([!]) can be written in an alternative way. 
Consider the U(N) currents 

J* = *7„A a * (4) 

with 

A° = hi (5) 
X j =t j , (6) 

P being the SU(N) generators normalized according to 



trtV = J<J y . (7) 

With these currents, one can define a non local, chiral invariant, U(N) Gross- 
Neveu model, with action given by 



S = J d 2 x*i0V- J <Px<PyJZ{x)V$(x-y)jfl(y), 

a,/? = 0,l,...,iV 2 -l (8) 

where V(^) are N 2 x A^ 2 symmetric matrices that weight the interaction. Taking 
N = 2, the actions (|l|) and (0) are equal provided that the matrices V( M ) can be 
written as 



V, 



O) 






V 

















(9) 



This non Abelian model was considered in Ref . ■ The bosonized effective 
action obtained by non Abelian bosonization led to a WZW functional, which 
was not easy to deal with in order to analyze the physical spectrum. In the 
next Sections we will show how this task is greatly simplified by starting from 
(P instead of (||) and combining path- integral Abelian bosonization and the self 
consistent harmonic approximation. 
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3 The equivalent bosonic action 

We start by considering the partition function 



Z=M J V^ a Vi> a e^ s , (10) 
where Af is a normalization constant. It is convenient to write 

S = So + S flip , (11) 

where 



.S'„= / d 2 x* a i#* a -^ [ <Px<Pyjftx)V$(x -y)j*(y) (12) 

and 



Smp = -9sj d 2 xd 2 y^ 1 (x) 7 ^ 2 (x)U {fl) (x-y)^ 2 (y) 7 ^ 1 (y). (13) 

The reason for this separation lies on the fact that So contains all the interaction 
terms that have separate chiral invariance for each species of fermions (spin 
states). They are Thirring-like interactions and will be treated in much the 
same way as the massless Thirring action, i.e. they will be transformed into 
free, non local action terms. Concerning the second term, it has no separate 
chiral invariance, and will be expanded in perturbative series like the mass term 
in the Thirring model. 

It is easy to show that the introduction of auxiliary vector fields allows 
to write 

Z=M' J V^> a V^ a VAl exp - J d 2 x § Q (0 + #t°) * Q - S[A] - S n 

where N' is a new constant (See || for details) and 

S[A] = \J d 2 xd 2 y (V^y" (x - y)A«( x )A^(y), (15) 
with (Vqj)^) defined through the following equation: 

d 2 V (v^y\x~y)V^(y-z) = 5^(x-z)6 ac . (16) 



(14) 
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We now decompose At in longitudinal and transverse pieces 



Al(x) = e^d u <t> a {x)+d^ a (x), (17) 

where 4> a and rj a are scalar fields. We also perform a change in the fermionic 
fields 



y a (x) = e-^^)-"^)]^) (18) 

whose Jacobian, as it is well-known, yields a kinetic term for the tfi a fields. One 
thus gets 

Z =N J Vx a T>x a V<f> a Vr] a e- Saa , (20) 
being S c g a sum of three pieces: 



where 



Sett — Sof + Sob + Sfiip (21) 



Sof = / d 2 x (x W + X W) , (22) 



9 - g2 
Sob - x 

z 



- J Sx [{d^f + (c^ 2 ) 2 ] + \ I d 2 xd 2 y (v { -') ab (x - y) 
x [e^ p d^ a (x)d p ^ b {y) + d^ a (x)d pV »(y) + 2e^M° \x)d^ b \y)} , (23) 



and Sfli p is the same SF interaction term already defined in equation (|l3|). 
Concerning this last term, from now on we shall restrict our study to the case 
of contact SF interactions: 



U(x - y) (0) - U(x - y) (1) = 5&(x - y), (24) 

and by performing a Fierz transformation followed by the chiral change defined 
in equations ( |l8| ) and (|l9|), we can write it in the form 



G 



Sflip = 2g s / d 2 x 



e -2 3 (0 1 ^)-l l + 75 x i . - 2 l-75 x 2 



• (25) 



Now we are ready to make an expansion of the partition function taking g s 
as perturbative parameter: 



= M J D0»2Ve-So B f^(_|jf)2 J (j[d 2 



i=i 



+ e 2 ^ 1 -^x 1 ^X 1 -X 2 ^X 2 



(26) 



OF 



where ( ) 0F means v.e.v. in a theory with action Sop. Only mean values involving 
an equal number of factors of the form 5X a (l + 7s)x a an d |X a (l — 75)x" ( n0 
sum over repeated indices is implied in these expressions) are non zero, and thus 
the partition function can be written in the form 



Z =J\f I V<f> a VT] a e- SoB 
(2 5s ) 2 " 



^ifgi_ J Jjj d 2 a ,. rf 2 y . e -2 s [^(x i )-0 2 (x i )-0 1 ( ! / i )+^(l/O] 



n=0 



\i=l 



1/ ^ 1 ~75. 



,i=l 



2 ~ v //v v " ' 2 

1 - 75 2^ x! + 75 



OF 
2/ 



x n^Y^^^v • (27) 

\*=1 ~ / OF 

The next step is to introduce two local massless scalar fields d a to be associ- 
ated with the free fermions yf an d X°- This trick allows to replace the fermionic 
mean values in the above equation by their ultraviolet regularized bosonic coun- 
terparts, which yields 
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Z = N j V(j) a Vri a e" S( " 

x£^|^ j /'jJ rf 2 a; . d 2 2/ . e -2 s [^(x i )-^(x i )-0 1 ( !/i )+^( ! / i )] N j 
n=0 ^ n '' \i=l ) 

x (— ) 4n /ne <VEF[ * 1(l4) "* 1( " <) A /jJ e ^K(x,)-^fe,)]\ _ (28) 

^ ' \»=1 / O,)? 1 \*=i / 0,# 2 

We now use the fact that the infrared divergences of the d a field's propagator 
provide a neutrality condition for v.e.v's of vertex operators. This allows to 
rearrange the perturbative series in a non-trivial way. One is then led to the 
completely bosonized action Sb os 

Sbos = Sob + J d 2 x ^(d^ a ) 2 - cos ^ig^ 1 - <f> 2 ) + v^tf 1 + ?? 2 )] J . 

(29) 

At this point it is interesting to observe that there is a set of changes of 
variables that allows to express the action in a very suggestive way. Indeed, 
writing 



-l^+tf 2 ) (30) 



9=±=(?-&) (31) 



i 1,2 = ^ c ±&) (32) 



V 1 - 2 = ^(nc±Vs) (33) 

where the plus(minus) sign corresponds to the pair (j) 1 , r\ x {<p 2 , if) , one sees that 
the field 6 completely decouples from the others and can then be integrated, 
and the bosonized action becomes 

Sbos = S c + Sg, (34) 

where 
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d 2 x (d^ c ) 2 + ljd 2 x d 2 y (Vfa) \x-y) 
x [^iiue lip d v (j) c {x)dp(j) c {y) + d^ c {x)d^ri c {y) + 2e llv d v (j) c {x)d l j,r] c (y)\ , (35) 



and 



x [e l±v e flp d u 4> s (x)dp^ s (y) + d tl rj s {x)d ll ri s {y) + 2e lll ,d v (j) s {x)d ll r] s {y)\ 



d 2 x 



9sA 2 



8tt9 



(36) 



with the functions ( V^* 



defined as 



2 y (v$) 1 (* - v)v^(jj -z) = sP\x z). 



(37) 



This, in turn, leads to the factorization of the partition function in the form 
Z = Z C Z S . This result is a clear manifestation of the so called spin-charge 
separation, a typical feature of Luttinger liquids |jl3f 0. Z c is the partition 
function associated to charge density excitations. It coincides with the spinless 
NLTM studied in reference Z s describes spin-density excitations. The action 
S s corresponds to a certain non local sine-Gordon model similar to the one 
previously considered in flio| |. In the next Section we shall derive an expression 
for the gap of its spectrum. 

Since our main goal is to analyze the spin-charge sector, from now on we 
shall focus our attention to S s . Going to momentum space (with the exception 
of the cosine term, whose transformed expression is not very illuminating) one 
has 



S, = 



d 2 p 

WT 2 



<j>{p)<j>(-p)A{p) + f)(p)f)(-p)B(p) + 4>(p)f)[-p)C(p) 
g s A 2 



J d 2 x cos (V8igcl> a + V&nO) (38) 



with 



vfa) 1 (p)pI + (v&) \p)pl 



(39) 
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B(P) = ~ 2 



(V^) 1 (p)p\ + (V ( s 0) ) \p)pl 



(40) 



C(p) 



PoPi- 



(41) 



In the above expressions <fi(p), flip) and #(p) are the Fourier transforms of 
<t> a {x), Vs(%) an d 6{x) respectively. 

4 The spectrum of the spin-charge sector 

In order to obtain the approximate spectrum for the spin-charge sector of 
our model, we will introduce a non local version of the self-consistent harmonic 
approximation |p"i"|| . Basically, this amounts to replacing the so called true action 
(p6|) by a trial action in which the cosine term is approximated as 



S A 2 



■ cos 



Q 2 

2 



8n6 



(42) 



where the parameter f2 of the trial action can be variationally determined (See 
the Appendix for details). Once this is done, it is straightforward to obtain 
the spectrum. Indeed, going to momentum space, and back to real frequencies, 
po = iu>, pi = k, the following equation is obtained: 



4tt0 2 

Our result for Q is 



fc 2 (l + V^(0))-^(l + V F (i)) 

TT w 7T w 



0. 



(43) 



O 2 = 



n A 2 



(44) 



where 



/o(fi) 



d 2 p 



^ 2 p 2 + £(p 2 V { \ )+ piVf 0) )+8n 



o 2 



(45) 



Equation (Q) is one of our main results. It gives, within the self-consistent 
harmonic approximation, a closed expression for the gap as functional of the 
potentials Vf^ip). 

At this stage it is interesting to consider the contact potential given by 
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V ( ' )(p) = V( a 1) (p) = l (46) 

In this case Io(tt) is infinite (See Appendix). After a suitable regularization 
one obtains 



A 2 7T 2 

This equation can be easily solved for A ^> ft and ^> A. The first case 
corresponds to <? s <C 1, and the result is 



(j_+_5_A0_ 
8tt^ 2 



(47) 



o 2 = ^4 



2 r 



85s 



Ws 2 



7r(l + <7 2 /7r) 



(48) 



Note that for the value g 2 = ir the gap posses a linear dependence on the 
coupling. It is interesting to observe that for this value, the action S s (see 
Appendix) corresponds to a free massive fermion, with mass proportional to g s . 
This is analog to the behavior found on the Luther-Emery line |l3) . Despite this 
similarity, the position of the line in our model differs slightly from theirs, but 
since the models are not exactly equal such a disagreement is not unexpected. 
However, we feel that the precise connection between both systems deserves a 
closer analysis that will be undertaken in a future work. 
The second case corresponds to g s 3> 1, and we obtain 



= (49) 

which is independent of g. This is expected because in this limit the contribution 
of the ordinary (without SF) terms are completely negligible. 

Let us stress that our gap equation (|44|) provides a new tool to check the 
validity of different potentials This is interesting because in condensed- 

matter applications the assumption of short-range electron-electron interactions 
works well for conductors in which the screening between adjacent chains reduces 
the range of interactions within one chain ]l6| . On the other hand, as the 
dimensionality of the system decreases charge screening effects are expected to 
become less important and the long-range interaction between electrons seems 
to play a central role in determining the properties of the system. This assertion 
seems to be confirmed by experiments in GaAs quantum wires Jl8| and quasi-ld 
conductors Jl7[ . We hope to report on numerical solutions of the gap equations 
for long-range potentials in a future contribution. 



11 



5 Summary 



In this work we have improved a non local version of the Thirring model 
which provides a tractable field-theoretical description of Luttinger liquids with 
spin. Indeed, in the context of non local Thirring-like theories previous treat- 
ments of SF interactions led to an involved non Abelian model (See for in- 
stance D). Specifically we have built an action based on two fermion species 
which allows to take into account SF interactions in an elegant and simpler 
way. Although our model was inspired by the one considered in Ref. p], it in- 
cludes interaction diagrams not contained in that previous work (the so called 
<74 diagrams in Solyom's terminology |T^| ). Besides, the theory we present has 
general bilocal potentials governing the interactions that do not flip spins. We 
parametrized this potentials in terms of the functions V c and V s which become 
associated to charge-density and spin-density dynamics respectively, once the 
spin-charge separation is made manifest after suitable changes of variables (See 
equations ( B2h and fl33|)). On the other hand, we could not get closed formu- 
lae for arbitrary distance-dependent SF potentials. Thus our analysis is only 
valid for local (contact) SF interactions. However, we were able to keep the 
general distance dependence of the ordinary (no SF) potentials up to the end 
of our computations. Under these conditions we got an effective bosonic action 
whose charge-density piece coincides with previously obtained descriptions of 
the forward-scattering spinless problem (p|. Concerning the more interesting 
spin-charge sector, we obtained, within the framework of the SCHA, a general 
equation which allows to have the gap in terms of arbitrary electron-electron 
potentials. Using this formula we could, for instance, estimate the effect of a 
long-range interaction on the gap. 
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6 Appendix: details of the SCHA method 

We shall give an sketch of the SCHA method. One usually starts from a 
partition function 



where Vfj, is a generic integration measure. An elementary manipulation leads 
to 




(50) 
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trial ( e 



for any trial action Stria! • Now, by means of the property 



.-(Strue-StrialV 



trial 

(51) 



( e -/)> e -</>, 



(52) 



for / real, and taking natural logarithm in equation (plf), we obtain Feynman's 
inequality Km 



111 -Ztruc > 111 Ztrial — (Struc — Stria! ) 



trial 



(53) 



In our case, the true action is the corresponding to the spin-charge dynamics, 
S s . However, before proceeding to the actual computation it is convenient to 
consider still another change of variables which diagonalizes the cuadratic part 
of S s . This change is given by 



A C 
t J *£! p 



g 2 p 2 

4B 2tt 



(54) 



= C iV^gCp 2 



4B 2tt 



(55) 



= -c 



and the resulting action, reads 



igp- 



(56) 



S a 



d 2 P 



C(pK(-p) 



Pi 



P&Vft+pjVfo ' T — r ' \v, 



(0) 



V, 



(1). 



■£(p)£(-p) 



P 2 + ^ (tiVfc +P\V( 0) ) I - 9 -^r jd 2 x cos(V8<) (57) 



where one sees that the £ and ^ fields become completely decoupled from £. 
Then it becomes apparent that the ^-dependent piece is the relevant one to 
consider in the present approximation. In other words, the relevant true and 
trial actions are given by 
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$true — 



(58) 



Atrial — 



d 2 p 



(59) 



where 



F(p)=p 2 + ^(p 2 V([ ) +p 2 



V (0) 



(60) 



The parameter fi can be determined by maximizing the right hand side of 
equation (p3) . In order to achieve this goal we first write 



In Ztnai = In / exp 



~ / d 2 xZ(x)(A0(x) 



= [ln(deti)~ 1/2 ] + const 



■trlnv± + const. 



(61) 
(62) 



where the operator A is defined, in Fourier space, by 



(AO(p) = [F(p) + 8irn 2 }t;(p). (63) 

It is then easy to get 

tr lni = V [ 4Ao HF(p) + 8vr0 2 ] (64) 
J [2nr 

where V is the volume (infinite) of the whole space J d 2 x. On the other hand, 
it is straightforward to compute (St rU c — Atrial)) by following, for instance, the 
steps explained in ref. pQ] . The result is 



(Strx 



tria l/ trial 



v 9sA 2 
= V — 5- exp 



-4tt 



V\-kVl 2 



d 2 p 1 
(2tt) 2 F(p) + 8irn 2 
d 2 p 



1 



(2tt) 2 F{p) + 8tt^ 2 



(65) 



Finally we can gather all the terms, and write them as 
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In Atrial ("Strue Atrial) trial 



v 



, exp (-47r/ (n)) + ^-h(V) - \h$l) 



where we have defined the integrals 



const. (66) 



d 2 p 
(2*)* 



\n[F(p) + Sirfl 2 



(67) 



with the formal properties 



d 2 p 



(68) 



(69) 



dJ-n(n) 



= -167r(n+ l)OJ_ n _i. 



(70) 



Now extrimizing expression ( pq) with respect to f2, and assuming that i_i(f2) is 
not zero (a condition that holds for most realistic potentials), we finally obtain 
the gap equation 



9/ 



e -47T/ (O) _ Q 



(71) 
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